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Abstract 



The two-dimensional Radon transform of the Wigner quasiprobability is introduced in 
canonical form and the functions playing a role in its inversion are discussed. The transfor- 
mation properties of this Radon transform with respect to displacement and squeezing of 
states are studied and it is shown that the last is equivalent to a symplectic transformation of 
the variables of the Radon transform with the contragredient matrix to the transformation 
of the variables in the Wigner quasiprobability. The reconstruction of the density operator 
from the Radon transform and the direct reconstruction of its Fock-state matrix elements 
and of its normally ordered moments are discussed. It is found that for finite-order moments 
the integration over the angle can be reduced to a finite sum over a discrete set of angles. 
The reconstruction of the Fock-state matrix elements from the normally ordered moments 
leads to a new representation of the pattern functions by convergent series over even or 
odd Hermite polynomials which is appropriate for practical calculations. The structure of 
the pattern functions as first derivatives of the products of normalizable and nonnormaliz- 
able eigenfunctions to the number operator is considered from the point of view of this new 
representation. 



1. Introduction 



The Radon transform of the Wigner quasiprobability is closely related to the expectation 
vàlues or densities formed with the eigenstates to the rotated canonical observables. They 
possess all properties of genuine probability densities and can be measured since recent time 
in quantum òptics of the radiation field by homodyne detection [1-4] . The field of problems 
of the reconstruction of the density operator from such or similar data is called quantum 
tomography. There arises the problem of the reconstruction of the density operator or of 
related quantities as matrix elements of the density operator or ordered moments from a 
complete set of tomographic data or of their approximate and, in some sense, optimal recon- 
struction from a more or less incomplete set of tomographic data. The complete tomographic 
reconstruction of the density operator from the Radon transform of the Wigner quasiprob- 
ability was prepared, as sometimes happens with things in history, in a wondrous manner 
by a theoretical paper of Vogel and Risken || shortly before the experimental realisation 
of the corresponding measurements in quantum òptics, although Vogel and Risken did not 
mention Radon transforms in this connection but speak about rotated quadrature phase. 
The transformation which is now called Radon transformation and which means the deter- 
mination of field functions from their integrated vàlues over some manifolds was introduced 
by Radon in 1917 but did not have a great impact on physics and also mathematics up to the 
fifties. Only, beginning from the sixties and prepared by mathematical papers as || and the 
requirements of different tomographic problems the Radon transformation and its inversion 
became interesting for a greater community that is reflected in monographs |7j and review 
articles [§. 

The theoretical development in quantum tomography in the last years went in the di- 
rection to determine more directly physical relevant parameters of the density operator as, 
for example, its matrix elements in the Fock-state basis [9-21] or in the position representa- 
tion p2| , [2"5| or to determine directly the normally ordered moments from tomographic data 
2\ p5|. Other sets of measurable quantities as the displaced Fock-state matrix elements 



of the density operator with fixed displacement parameter and varied excitation numbers 



are investigated for its appropriateness to the reconstruction of the density operator [26 



Another direction is to study the influence of imperfect measurements onto the Radon trans- 
form of the Wigner quasiprobability. If this influence consists only in a smoothing of the 
Wigner quasiprobability and of its Radon transform resulting in its convolution with a Gaus- 
sian function then there is, in principle, no loss of information and one has to take one of 
the smoothed quasiprobabilities instead of the Wigner quasiprobability. This was taken into 
account in some of the first papers about reconstruction ( e.g., |TJ[] ) and was prepared 
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already by the paper of Vogel and Risken |J starting from the rotated quadrature compo- 
nents of the s-orderd quasiprobabilities. Another problem arises when there is a real loss 
or incomplete knowledge of information. Then the reconstruction is no more unique and 
one needs a principle for an optimal reconstruction under such conditions. This problem 
was recently studied from a general point of view in [27J where the incomplete knowledge 
was called observation level and this could become important in quantum tomography in 
the next time. The necessary principle is Jaynes' principle of maximal entropy and the best 
choice for the entropy in quantum mechanics is, in our opinion, for some formal reasons the 
Von-Neumann entropy. 

The direct reconstruction of the matrix elements of the density operator in the Fock- 
state basis initiated by D'Ariano, Macchiavello and Paris || leads to the integration over 
the angle and over the position of the rotated quadrature components multiplied by specific 
functions of the angle and of the position coordinate. These auxiliary functions split into 
simple angle-dependent phase functions and more complicated position-dependent functions 
which were called pattern functions [12-14] and were calculated in a representat ion by the 
functions of the parabòlic cylinder. Then it was found that these pattern functions can 
be represented as first derivatives of products of the normalizable eigenfunction with one 
nonnormalizable eigenfunction of the number operator to eigenvalues corresponding to the 
numbers in the considered Fock-state matrix elements |fL8|| . The best way to recognize this 



special structure of the pattern functions is to consider the differential equations for the 
products of the eigenfunctions ( Hermite functions for the harmònic oscillator ) and to 
establish the corresponding orthogonality relations by Standard methods of the theory of 
ordinary differential equations that means with the help of the adjoint equations ]T9| , [2üj. 

The usual restriction of the representation of the Radon transform in quantum òptics 
to the dependence on two variables ( rotation angle and line coordinate ) is unfavourable 
when considering the transformation properties with respect to squeezing of the states which 
results in symplectic transformat ions of the arguments of the Wigner quasiprobability and 



of its Radon transform. Therefore, in |25| was chosen a more general concept for the Radon 
transform as the starting point and it was explicitly calculated this Radon transform for 
squeezed coherent states. The same more general concept was taken in [[28|, ^9] and called 
symplectic tomography and since the Radon transform contains complete information about 
the quantum-mechanical state the bàsic equations of quantum mechanics were reformulated 
for the Radon transform. This has the advantage that one works directly with measurable 
quantities which can be considered as genuine probability densities, at least, for each separate 
choice of the angle. 

In the present paper we introduce in section 2 the two-dimensional Radon transform, as it 
seems to us, in its most rational form which we call canonical form and consider its inversion 
where our intention is to clarify the meaning of some mathematical aspects which play a 
role in the inversion of two-dimensional Radon transforms and to give a new representations 
of this inversion. In section 3 we consider the transformation properties of the Wigner 
quasiprobability and of its Radon and Fourier transforms under displacement and squeezing 
transformations of initial states. The developed formulae, although highly technical and 
partially complicated, seem to us as very important for practical applications and theoretical 
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calculations of the influence of squeezing that will be demonstrated for squeezed coherent 
states. Section 4 is devoted to the introduction and calculation of the pattern functions [12- 
14] for the direct reconstruction of the Fock-state matrix elements from tomographic data 
and follows with weak deviations the main stream of considerations in the literature about 
quantum tomography but this section is also important as a preparation of the content of 
the next two sections. In section 5 we consider in a short form the direct reconstruction of 
the normally ordered moments from complete tomographic data without integration over the 
angle as recently published in our paper []25] and give some new formulae not represented 
there. Starting from the reconstruction of the density operator via the normally ordered 
moments we derive in section 6 an essentially new representation of the pattern functions 
for the reconstruction of the Fock-state matrix elements of the density operator in form 
of series over Hermite polynomials of even or odd order which is an alternative for the 
practical calculation of these functions to the existing representat ions by functions of the 
parabòlic cylinder or by first derivatives of the product of the normalizable with a certain 
nonnormalizable eigenfunctions of the number operator in "position" representation. We 
find in this section a nonuniqueness of the pattern functions for the reconstruction of the 
Fock-state matrix elements. In section 7 we derive the specific structure of the pattern 
functions as derivatives of a product of a normalizable and a nonnormalizable wave function 
from the fourth-order differential equation for products of Hermite functions and its adjoint 
differential equation. Here we come again across with the mentioned nonuniqueness of the 
pattern functions. 
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2. Two-dimensional Radon and 
Fourier transforms and their inversion 



The Wigner quasiprobability is the best compromise in quantum mechanics for a phase-space 
description of a quantum-mechanical state in analogy to the phase-space description in clas- 
sical mechanics and statistics [30-38]. We consider one field mode corresponding to one 
mechanical degree of freedom and use a pair of canonical coordinates (q, p) or complex coor- 
dinates (a, a*) and corresponding pairs of canonical operators (Q, P) or boson annihilation 
and creation operator (a, at) in the following way 

q + ip t q-ip Q + iP t Q-iP i * dq A dp 
a = — a = — a— — a 1 = — -da A da = — . (2.1) 

V2h v 7 ^ V2h 2 2h K J 

In this section we consider the Radon transform of the Wigner quasiprobability its connection 
to the Fourier transform of the Wigner quasiprobability and the inversion of these transfor- 
mations. All considerations of this section are not specific for the Wigner quasiprobability 
and one can exchange the symbol for the Wigner quasiprobability by the symbol for an arbi- 
trary other function of two variables over a plane, for example, another quasiprobability and 
one obtains in this way the general properties of two-dimensional Radon transformations. 
Specific properties of the Wigner quasiprobability and its Radon and Fourier transforms are 
considered from section 3 on. We use here the representation by the real variables {q,p). All 
these relations can be easily translated into a representation by complex variables. 

The Radon transform W(u,v;c) of the Wigner quasiprobability W(q,p) can be defined 
in the following canonical form ( compare [6-8] ) 

W(u,v;c) = j dq A dp5(c- uq - vp)W(q,p). (2.2) 

The normalization of the Radon transform is closely connected to the normalization of the 
Wigner quasiprobability in the following way 

/ dcW(u,v;c) = / dqAdpW(q,p) = 1. (2.3) 

J —oo J 

Due to the relation 

1 . 

W(fxu, fJ,v;(AC) — |— tW(u, v;c), fj, e R, (2.4) 
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where /i is an arbitrary real number this Radon transform depends effectively only on two 
continuous variables. The presence of the delta function under the two-dimensional integral 
in (2.2) restricts the integrations to the one-dimensional objects 

uq + vp = c. (2.5) 

These are the equations for straight lines with (u, v) as a normal vector to the lines and 
c/ \/u 2 + v 2 as a measure for the orthogonal ( nearest ) distance of the line to the coordinate 
origin. However, this "oriented distance" can take on all real vàlues contrary to positively def- 
inite distances. The coordinates (u,v; c) are homogeneous line coordinates in the dual plane 
to the (ç, p) -plane. In case of the gener alization of the Radon transform to iV-dimensional 
spaces the corresponding (N — l)-dimensional objects over which the field is integrated are 
the (N — l)-dimensional hyperplanes but this will be not considered here. Due to Eq.(2.2) 
the full information of the Radon transform is already contained in W^(cosy3, sin<^; q v ) where 
(u, v) = (cos <f, sin <f) is now the normal unit vector to the straight lines. This reduced Radon 
transform is mostly denoted by w(ip, q v ) but we will see in the next section that such a nota- 
tion, although shorter, is unfavourable for the discussion of the transformation properties of 
the Radon transform. We call the Radon transform as introduced in Eq.(2.2) the canonical 
form. 

The Radon transform W(u, v;c) is closely related to the Fourier transform W(u,v) of 
the Wigner quasiprobability W(q,p) which can be defined by 

W(u,v) = j dq A dp exp{—i(uq + vp)}W(q,p), (2.6) 

with the inversion 

W(q,p) = -j—^ J duAdvexp{i(uq + vp)}W(u,v), (2.7) 
and with the normalization 

W(0,0) = J dqAdpW(q,p) = 1. (2.8) 
The relation of the Radon transform to the Fourier transform is given by 

W(u,v) = J dcexp(-ibc)W í-,-;c) , b e R, (2.9) 
with arbitrary real numbers b and its inversion is given by 

1 r+oo 

W(u,v;c) = — db exp{ibc)W{bu,bv). (2.10) 

The full inversion of the two-dimensional Radon transform can be made in two steps. The 
first step is the transition from the Radon transform to the Fourier transform and the second 
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step is the inversion of the Fourier transform. If we introduce polar coordinates (r, ip) instead 
of (u,v) the following integral can be accomplished 

/+oo Q / r+oo rO 

dr \r\ exp (—ixr) = i— ( / dr exp (—ixr) — dr exp (—ixr) 
-oo OX \Jo J-oo 



d 



+ 



dx \i(x — iO) i(x + iO) J 
= ïS-V-, (2-11) 

OX X 

where V in front of 1/x means that integrals over the singularity of this function must be 
taken in the sense of Cauchy's principal value. Then we find 

r 2 r n r+oo (d i \ " 

W(q,p) = -^ dtp dc \^-' p 7 ; : r )W(r cos (p,r sm ep; c), 

2tt jo J-oo \oc c — r (q cos tp + p sm p) J 

(2.12) 

with arbitrary real fixed r§. The integral over c can be represented after partial integration 
as the Hilbert transform of the derivative of the Radon transform with respect to the variable 
c. A more symmetrical way of the representation of the inversion formula can be obtained if 
first the integration over c is made in the two-step inversion via the Fourier transform. One 
finds 

1 f ( 1 - 

W(q,p) = — du Adv [V \c\W(u,v;c), c G R, (2.13) 

2n 2 J \ c — uq — vp J 

where the upper index (1) denotes the first derivative and where c is an arbitrary real number. 

It seems that the meaning of (d/dx)V(l/x) is not very clear in the literature about 
quantum tomography. In particular, one cannot write — V (1/x 2 ) as the result of this differ- 
entiation. Let us therefore discuss this more in detail. One has to consider (d/dx)V(l/x) as 
a generalized function or linear functional and has to define all its derivatives as generalized 
derivatives in the sense of the theory of generalized functions. We introduce here the symbol 
1Z for the canonical regularization of a singular function ( Gel'fand and Shilov JÏÏJ use the 



symbol k.r. in the Russian and CR in the English edition for canonical regularization ). 
Canonical regularizations of different singular functions are synchronized in the sense that 
with respect to additions of functions, multiplication of functions by numbers and by well- 
behaved smooth functions and differentiations the regularization symbol can be changed in 
its order with these operations. The canonical regularization of the singular function 1/x 
is identical with the principal value V(l/x). If <p(x) denotes arbitrary basis functions then 
one can define 71(1 /x) by the following linear functional ( we do not specify in this physical 
context the space of basis functions which must be sufficiently smooth, in particular, at 
x = and sufficiently decreasing in infinity together with their derivatives ) 

n 1 -, p(x)) = r d x (n 1 -) p(x) = f-js 

X ) J-oo V xj Jo X 

í +CO dx (v- )<p(x) = lim / dx^-, TZ- = V- = ^-\og\x\. (2.14) 

J-oo V XJ z^J\x\>e X X X OX 
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The function log |x| is a regular function because it is locally integrable at x = and all 
canonical regularizations of power functions l/x n can be defined as derivatives of log |x| 
( Gel'fand et al. || omit the regularization symbol TZ when discussing the functions 
7Z(l/x n ) but understand these functions in the sense of canonical regularization ). The 
canonical regularization of 1/x 2 can be defined by 



X 2 ' / J~oo V X 2 J JO X 2 

1 dl d 2 

n- = -j-n- = -— log x . (2.15) 

X A OX X ox z 

The proof of this formula can be given by the following chain of transformations including 
partial integration and using (2.14) 

k\,<p(x)) = -(-n-,^\x) 

x z / V x 







y(g) + <p(-x) - 2<p(0) _ f (p(x) + <p(-x) - 2<p(0) \ 



dx 



—x 2 x 



) x=0 

(2.16) 



The last expression in braces taken at x = vanishes and thus Eq.(2.15) is proved. All other 
proofs of canonical regularizations can be made in analogy to this proof which is here given 
as an example for such kind of proofs. The meaning of {V{l/x)) (l) in the inversion formulas 
for the Radon transform is therefore that we have to substitute it by — 71(1/ x 2 ) considered 
in the sense of Eq.(2.15). One may think that this does not give the right sign but it is not 
so. As an example we consider (p(x) = exp(— x 2 ) 

1 2 \ [+°° 2(exp(-x 2 ) - 1) 



1Z—, exp(— x 2 )] = / dx 
x l ) Jo 



X 2 



- 3C : 



dx 4X6XP( ^ = (li-, -2xexp(-x 2 )) = -2v^F. (2.17) 

o x V x ) 

It could seem that lZ{l/x 2 ) and exp (— x 2 ) are both positively definite functions ( exp (— x 2 ) 
is really positively definite ) but, nevertheless, the value of the linear functional in (2.17) for 
the considered basis function gives a negative number as shown. This means that 71(1 /x 2 ) 
is not a positively definite function. 

The two-dimensional Radon transformation (2.2) together with its inversion (2.13) is 
equivalent to the following resolution of the two-dimensional identity operator 

|r ' r duAdv(ll-, j 5{c-uq'- vp') = ó(q - q')5(p - p'). (2.18) 



27r 2 J \ (c — uq — vp)' 
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All the seeming difíiculties with the inverse two-dimensional Radon transform can be avoided 
if the integration leading to the inversion in two steps via the Fourier transform as the 
intermediate step is not made at an early stage of the calculations. 
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3. Transformation properties of the 
Radon transform of the Wigner 
quasiprobability 



The Wigner quasiprobability can be obtained by forming the trace over the product of the 
density operator with a transition operator. Among all quasiprobabilities it has the unique 
property that the operators for the transition from the density operator to the quasiprob- 
ability and back from the quasiprobability to the density operator via integration are the 
same. The transition operators possess many equivalent representations. For our purposes 
we use the following convenient representation (e.g., 



W(q,p) = (gexp {^Q^~ q ~ ) <W(p), J dq A dpW{q,p) = 1, 

Q = 2Tm J dqAdpW(q,p)exp(^-Q-^-- P-^j 5(q)5(p), (A) = Trace{A}, (3.1) 

where we use the very rational notation (A) for the trace of an operator A. We now inves- 
tigate the transformation properties of the Radon transform of the Wigner quasiprobability 
under displacement and squeezing transformations of the density operator. 

Suppose that there is given a density operator g corresponding to a Wigner quasiprob- 
ability Wo(q,p)- Let us make a displacement of the state described by the density operator 
Qq by using the displacement operator D(q,p) in the following way 

g = D(q,p)g (D(q,p))\ D(q,p) = exp [~{qP ~pQ)} = (D{-q, -p)) ] . (3.2) 
Then we find for the Wigner quasiprobability W(q,p) of the displaced state g 
W(q,p) = ^g (D(q,p)yexp^-Q^-P^D(q,p)y(q)6(p) 

= (qo exp j-(Q + ql) ^ - (P + Pl)^j S(q)5(p) = W (q -q,p- p), (3.3) 

meaning that the new Wigner quasiprobability is simply the displaced old Wigner quasiprob- 
ability in the phase plane. This property is also true for all other quasiprobabilities after 



9 



displacement of the state due to the displacement structure of the quasiprobabilities them- 
selves ||38|| . The Radon transform of the primary and of the transformed state are then 
connected by a displacement of the variable c 



W(u, v; c) = Wo(u, v;c — uq — vp), 
whereas the Fourier transforms are connected by 

W(u, v) = exp{—i(uq + vp)}W Q (u, v) 



(3.4) 



(3.5) 



that means by multiplication with a phase factor depending on (u, v) and on the displacement 

(q,P)- 

Our next aim is to consider the influence of unitary squeezing transformations of a state 
described by the density operator onto the Wigner quasiprobability and its Radon and 
Fourier transforms. We prepare this by considering first some bàsic formulae for the, in 
general, nonunitary squeezing operator S(Ç, rj, () which is defined by 



S(Ç,V>Q = exp||a 2 + i|(aa + + a + a) - ^a t2 j 



exp 



% 

2h 



^-0 + 2 V + í±í (QP + PQ) + ^-^ + 2 V 



(3.6) 



The transformation of the canonical operators (Q, P) by the general squeezing operator is 
given by ( in representation by the boson operators (a, a^) this is given in [^] and more 
specially in ) 



S(t,Ti,Q(Q,P)(S(t,Ti,CT l 
(Q,P) 



( h ^ ~ C) ~ 271 sh£ ^ 

° £ 2 e ' 2 e 



\ 



iíjj -Q + 277 she g + Ç she 

2 e ' C £ 2 e 



= (Q,P) 



a, (3 \ 


a, (3 


7, à )> 


7, 5 



(3.7) 



with the inverse correspondence 
a — i/3 — i'y — 5 Arshii 



/3 — 7 Arsh# 
2 ïT~ ' 



a + + 27 — <5 Arsh# 
2 ~ ' 



\ 



a + <T 



1 = ±she. 



(3.8) 



This is the fundamental two-dimensional representation of the SL(2, C) ~ 5*^(2, C)-group 
which in case of unitary squeezing operators £>(£'*, r/ = r)'*,Ç') reduces two the SL(2,R) ~ 
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Sp(2, R) ~ SU(1, l)-group |4^, |13| . The matrices of this representat ion in (3.7) are com- 
plex unimodular or symplectic but in general not unitary. Symplectic transformations in 
spaces of even dimensionality are linear canonical transformations. They preserve an anti- 
symmetric 2-form which can be nondegenerate only in spaces of even dimensionality whereas 
unimodular transformations preserve a fully antisymmetric n-form ( volume form ) in ar- 
bitrary n-dimensional spaces. Only in two-dimensional spaces one has coincidence between 
symplectic and unimodular transformations. The notions of symplectic transformations and 
symplectic groups were introduced by Weyl JBJ] . We mention here that all problems concern- 



ing the multiplication of squeezing operators or the decomposition of squeezing operators 
into products with certain ordering are solvable by applying the fundamental representation 
of the two-dimensional symplectic or unimodular group in Eqs.(3.7) and (3.8) ( see also f|5| 
and H ). 

In the special case 

f = C = 0, V = V* = <f, S(0, ip, 0) = exp U~ \ R(ip), R(ip) = exp [iipa^aj , (3.9) 
this operation reduces to rotations according to 

R( ( p)(Q,P)(R( í fi)y = (Q,P) ( ™% ) = {QcoB<p + Pmn<p,-Qfàn<p + PaoB<p), 

(3.10) 

where R(<p) denotes the rotation operator. 

We now consider the following squeezing transformation of the primary state g by a 
unitary squeezing operator £>(£'*, 7/ = r]'*,Ç') according to 

q = S(C, rf, Ç')Qo(S(Ç'*,r]', C))\ (v' = v'*), (3.11) 

where the unitary squeezing operator S(Ç'*,r)' = rj'*,Ç') is a special case of the, in gen- 
eral, nonunitary squeezing operator S(Ç,r),() defined in Eq.(3.6). By using the trans- 
formation of the bàsic canonical operators (Q,P) given in (3.7) and taking into account 
(S(C*, i = i*, C')) 1 = S(-C'*, -i = -rf*, -C) one finds 

W(q,p) = /^ o (^(C'^r 7 ^C0) t expf-g^-P^^(C / ^r 7 ^C0\%)í(p) 
= L exp L(Q6 - P 7 )^ - (-QP + Pa)-^V\8(q)5(p) 



[^(-Q^-P^))s {qW )) 



(q'=aq+yp, p'=(3q+5p) 

= W (aq + W ,f3q + 5p). (3.12) 

This means that the Wigner quasiprobability W(q,p) for the state after squeezing can be 
obtained from the initial Wigner quasiprobability Wo(q,p) by an argument transformation 
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with a real unimodular or symplectic matrix in the following way 



W(q, p) = W (aq + jp, f3q + Sp) = W (ç', p') , (</, p) = (q, p) 



a, j3 
7, 5 



7, ó 



-Ç*)-2r)' she' \ 



\ 



C'C* - v 



r2 



2 e' 1 
-i(C - C'*) + 2r/' she' 

rf = V). 



che' 



2 e' 



(3.13) 



This result was originally obtained by Ekert and Knight g6|, 47]. 

The Radon transform W(u, v;c) of the Wigner quasiprobability after squeezing of the 
state can be obtained from the Radon transform W (u, v\ c) of the Wigner quasiprobability 
of the initial state by ( the double meaning of the symbol 5 in the following formula can 
hardly lead to confusion ) 



W(u,v;c) = I dq A dpó(c — uq — vp)W(q,p) 

dq A dp5(c — uq — vp)Wo(aq + 7p, (3q + 5p) 

dq' A dp' ó(c - (5u - (3v)q' - (—ju + av)p'^W (q' ', p) 

= Wq (5u — /3v, — r yu + av; c) , 

which can be also written in the form 

W(u, v; c) = Wç> (5u — fiv, —^u + av; c) = Wo(u', v'; c), 

-i / \ 



(3.14) 



u 
v 1 



S, 
"7, 



a 



u 
v 



a, P 
7, 5 



u 

V 



(3.15) 



The transformation of the variables (u, v) in the Radon transform is made with the contra- 
gredient ( inverse transposed ) matrix to the matrix of the transformation of the canonical 
variables (q,p) if we consider (u, v) as a row vector and by the inverse matrix if we consider 
(u, v) as a column vector. The same is true for the transformation of the variables (u, v) of 
the Fourier transform for which one now easily finds 



W(u, v) = W (5u - Pv, -ju + av) = W (u', v'). 



(3.16) 



It could appear that it is more simple to derive first this result for the Fourier transform of the 
Wigner quasiprobability and then the corresponding result for the Wigner quasiprobability 
itself but, really, both derivations are of the same degree of difficulty and very similar. 

The combination uq + vp is invariant with respect to the considered unitary squeezing 
transformations of states 



ii, ii 
u q + v p 



(q',p') 



u' 
v' 



u 

V 



uq + vp, 



(3.17) 
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but not with respect to displacements of states. One can make both displacement and 
squeezing transformat ions of the initial state but then the result for the transformed Wigner 
quasiprobability and its Radon and Fourier transforms depends on the order of these op- 
erations. If we make first the squeezing operation and after this the displacement then we 
can immediately apply the derived transformation properties. However, if we make first the 
displacement and after this the squeezing of the state then we have to take into account 
that the displacement parameters must be transformed by the real unimodular or symplec- 
tic matrices in the same way as the corresponding coordinates. We do not explicitly write 
this down. Usually, the Radon transform W(u, v; c) of the Wigner quasiprobability W(q,p) 
is only considered in the special form W(cos(p, simp; q v ) and denoted by w((p,q v ) but it is 
difficult to express such a fundamental transformation property as in Eq.(3.15) in the last 
notation. Let us also mention that the transformation properties of all other quasiproba- 
bilities and their Radon and Fourier transforms with respect to the considered squeezing 
transformations of states are not so simple as for the Wigner quasiprobability and we will 
not consider this here. 

Mancini, Man'ko and Tombesi [ 2"8"| , ^ consider practically the same object W(u, v;c) 
which is the Radon transform of the Wigner quasiprobability ( slightly extended by a further 
displacement ) and call this symplectic tomography. Since W(u, v\ c) contains the complete 
information about the state one can transform the equations of time evolution for the state 
from any other quantity with complete information, for example, for the density operator 
or for the Wigner quasiprobability to an equation for the Radon transform. A certain 
advantage of this approach is that we consider then directly the measurable quantities which 
are positively definite ( or semidefinite ) and can be considered for every fixed (u, v) as 
genuine probability densities |28, |29|]. The noncommutativity of the bàsic observables in 



quantum mechanics appears then as the, in general, noncompatibility of all these probability 
densities for essentially different (u, v) on a classical level that means as marginals or Radon 
transforms of purely classical distribution functions f(q,p). 

Let us shortly consider the transformation of the Wigner quasiprobability in the complex 
representation with respect to unitary squeezing transformations. The corresponding formula 
to (3.7) is then 

S(t,Tí,Q{a,ai){S{t,ii,0)- 1 

she 

É — 



a, a 



^ che — ir] , 



she 



she 



V 



C , che + irj / 

e e ' 



with the correspondences 




k, X 
/i, v 



K, X 
/i, V 



tt-v 2 



(3.18) 



/ a + ij3 — + 5 a 
\ 



i/3 — ry — 6 \ 

2 ' 2 

a + i(3 + ry — 5 a — i(3 + ry + 5 

"9 J 



k + v = a + 5, 
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a, j3 
7, 5 



/ K + X + fM + V .K — À + /i — V \ 

2 ' ~ % 2 

.K + X—fl — U K — X — 11 + V 



\ 



(3.19) 



/ 



2 1 2 

Then one finds the transformed Wigner quasiprobability after the unitary squeezing trans- 
formation of the state according to (3.11) in analogy to (3.12) 

W(a, a*) = W (v*a + tia* , ffa + vot) = W (a', oi"), (3.20) 

with the following specified unimodular or symplectic matrix due to the unitarity of the 
considered state transformation 

she' 



/X, V 



i 



\ 



che' — ir] 



, she' 



c 



, she' 



C 

f 

che f + ir]' 



she' 



) 



(3-21) 

The propagation of the transformation properties (3.20) of the Wigner quasiprobability in 
complex representat ion to its Radon and Fourier transforms in complex representation is 
easily to find if one introduces a certain convention for the defmition of these complex 
representations. We will not do this here. 

Let us now consider the more special unitary squeezing operator S(('*,0,(') as it is 
mostly used in the theory of squeezing [48-58]. According to (3.6) we have 



S(Ç'*,0,£) 



JC* 2 C t a' 
exp< — a - -a' 



cxp 



C th|Cl 



,t2 



2 K'l 



l H^rw-iN+L·) ( ('* th|C'| 

(ch|C |) 2 ; exp I — — — 



ICP 



* C 2 

exp — a 



N = a f a, (3.22) 



where there is additionally given an almost normally ordered decomposition into factors ( 
the normal ordering of the central factor is mostly not necessary in applications but could 
be easily made ). Furthermore, it is made a new parametrization of the squeezing by the 
complex variable ( characterized by 



C = C'^, C' = C^p, ICI = th|C'|, |C'l = Arth|c|. 



(3.23) 



The complex parameter ( describes the squeezing operation in the nonunitary approach 
pï| , |59[| For convenience of applications of the preceding formulae we give here the explicit 
form of the corresponding specified matrices 

C' + C* sh|C1 
a, P \ _ -M-*-— 2 icf' 



7, 



— i 



X'-C* Bh|C'| \ 

í 



V 



X'-C* sh|ç'| 
1 2 K'l 



ch|C'| 



2 K'l 
C' + C* sh|c'| 
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i-ICP 



V 



2 

.c-c* 



. c-c \ 

2 

c+c 



I/' 



( chici, 

, shiC'l 



K'l 



ch|C| 



i, c 



} 



^i-ICi 2 vc, i 



(3.24) 



Now we can apply all derived transformation formulae for the Wigner quasiprobability and 
its Radon and Fourier transforms and have to specify the corresponding matrix elements by 
the expressions given here. 

The substitution of the complex squeezing parameter Ç' by Ç according to Eq.(3.22) is 
not only a formal substitution but has a deep physical content. If one acts with squeezing 
and displacement operators onto pure states \ip) one can represent \ip) by its Fock-state 
representation. Then one can represent the Fock states by the action of powers of the 
creation operator onto the vacuum state |0) that means by an excitation operator. Then 
one can bring the squeezing operator of the form (3.22) by application of the commutation 
rules with the excitation operator in the position that it acts in unchanged form directly 
onto the vacuum state. By using the ordered decomposition given in (3.22) the operator 

exp (((* /2)a ) acts onto the vacuum state as the unity operator whereas (vi — |C| 2 ) 

multiplies it with a normalization factor (1 — ICI 2 ) 1 ^ 4 - The essential action of the squeezing 
operator onto the vacuum state can be then substituted by the action of the nonunitary 
squeezing operator exp (— ((/2)a^ 2 ^j. After this the changed excitation and displacement 
operators act onto the obtained state that must be considered for the concrete case. This 
procedure should be also possible for mixed states or genuine density operators. Hence the 
squeezing parameter ( in the nonunitary approch possesses a concrete physical meaning. In 
case of the more general unitary squeezing operator S(Ç'*,rf = rf*,Ç') one has 1 complex 
parameter (' and 1 real parameter rf that is more than the only 1 complex parameter ( in the 
nonunitary approach and there remains a certain freedom in the choice of the correspondence 

We do not discuss this here. 



between unitary and nonunitary approaches [40 



As an example we consider the vacuum state qq ■ 
and its Radon and Fourier transforms are 



|0)(0|. The Wigner quasiprobability 



W (u,v;c) 



Wq(u, v) = exp 



exp 



Hir(u 2 + v 2 ) 
Ti(u 2 + v 2 ) 



fi(u 2 + V 2 ) 



(3.25) 



After squeezing of the state with the special unitary squeezing operator in Eq.(3.23) and 
succeeding displacement one obtains squeezed coherent states and by using Eqs.(3.3)-(3.5) 
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and (3.13)-(3.16) as well as (3.24) one finds the following Wigner quasiprobability and its 
Radon and Fourier transforms |f25|| 



w , v 1 [ |1 + C| 2 (9 - q) 2 + |1 - C\ 2 (P - P? - i(C - C)2(q ~ q)(p-p) 
W(u,v; c) 



ICP 



\ h-n (|i - c\ 2 u 2 + |i + Cl 2 *; 2 + i(C - C*)2w) 

f (1 - ICI 2 ) (c-uq-vp) 2 



x exp 



& (|i - Cl 2 « 2 + |i + Cl 2 « 2 + i(C - C*)2w) J ' 

W^(ií, u) = exp {— i{uq + vp)} 

^^MIi-CIV + li + CIV + iK-C-)^)!, (326) 

We mention the following complex factorizations of the numerators and denominators in 
the exponents which are possible with rational coeficients in ( and (* only for the Wigner 
quasiprobability and its Radon and Fourier transforms but not for other important quasiprob- 
abilities of these states such as, for example, the coherent-state quasiprobability. 

|i + Cl V + |i - Cl V - i(C - C)2q P = {(i + C)q + »(i - CM{(i + C)q - i(i - CM, 
|i - Cl V + |i + Cl V + z(C - C)2uv = {(i - ()u + i(i + C)«}{(1 - C> - <(i + C)«}- 

(3.27) 

They express a certain duality of the transformation of the variables in the Wigner quasiprob- 
ability on one side and its Radon and Fourier transform on the other side. 

For convenience we give in the following formula the connection between our parameters 
q, p and C for squeezed coherent states and the expectation vàlues of the canonical operators 
(Q, P) and their dispersions and symmetrical correlation ( A = (qA), AQ = Q — QI, AP = 
P-PI ) 

^"Cl 2 7X7^ ^ll + CI 2 



Q = q, P = p, (AQ) 2 = -i (AP) 



21- |C| 2 ' 2 1 - |C| 2 ' 



.n C-C 



-(AQAP + APAQ ) = <2 rq^p- ( 3 - 28 ) 

These 5 real parameters can be used for the unique characterization of squeezed coherent 
states. Because of 



i _______ _ __ v 2 h 



2 



(AQ(<p)Y (AP(<p)y--(AQ(<p)AP(tp) + AP(<p)AQ(tp)) = — , (3.29) 

for arbitrary ip they effectively reduce to 4 independent parameters for squeezed coherent 
states. However, for general Gaussian Wigner quasiprobabilities corresponding to displaced 
and squeezed thermal states these parameters become independent and (3.39) has to be 
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substituted by an inequality corresponding to a modification of the usual uncertainty rela- 
tions by inclusion of the uncertainty correlation [60-62]. The maximal and minimal vàlues 



of y (AQ((p)) for varying ip are determined by 



(AQ( Vmax )y 

exp(i2Lp max ) = - 



fi 1 + lCI 

N 2 1- icr 

exp(i2(p min ) 



c 



+ 



icr 



fi í-ICI 

>|2 1 + ICI' 



(3.30) 



These uncertainties can be considered as one of the possible definitions of the widths of the 
principal axes of the squeezing ellipse. For the corresponding angles ip ext determined by 
(3.30) the symmetrical correlation vanishes, i.e. 



AQ(<p ext )AP(<p ext ) + AP(<p ext )AQ(<p ext ) = 0. 



(3.31) 



In an analogous manner one can consider squeezing and displacement of other initial 
states, for example, Fock states, thermal states or Schròdinger cat states. 
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4. Reconstruction of the density 
operator from the Radon transform 



In this section we consider the reconstruction of the density operator and of its matrix 
elements in the Fock-state basis from the Radon transform of the Wigner quasiprobability. 
Let us first establish the following relation between the Radon transform and the rotated 
marginals of the Wigner quasiprobability 

W(cosip, smip; q) = (q\(R(ip)) J gR(tp)\q) = (q; ip\g\q; <p), (4.1) 

where \q; <p) denotes the orthonormalized eigenstates of the rotated canonical operator Q(<p) 
to eigenvalues q in the following way ( see Eq.(3.10) ) 

Q(<p)\q; tp) = R((p)Q(R((p))*R((p)\q) = q\q; <p), Q\q) = q\q), \q; <p) = R(<p)\q), 

/+oo 
dq\q;(p)(q;<p\ = I. (4.2) 
-oo 

To prové Eq.(4.1) we first immediately find the special result 

d P W(q,p) = W(l,0;q). (4.3) 

-oo 

Then, we see from Eq.(4.1) that (q; <p\g\q; (f) corresponds to the Wigner quasiprobability of a 
rotated density operator about an angle — cp that means to R(—(p)g(R(—ip)) !í and according 
to Eqs.(3.10), (3.12), (3.14) or the definition of the Radon transform in Eq.(2.2) one obtains 

/+oo 
dp W^gcos (p — psin ip, q sin ip + pcos tp) = ^(cos tp, sin ip; q). (4.4) 
-oo 

Thus the relation (4.1) is proved. From the well-known position representation (q\n) of the 
Fock states |n) it follows 



With the summation formula of Mehler í|63 



- z n 1 j 2xyz- (x 2 + y 2 )z 2 
> :H n {x)H n {y) = —. exp < 



z 2 



(4.6) 
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which is not difficult to prové [ÏO] one calculates by Fock-state expansion the more general 



scalar product ( Green's function to the wave equation for the harmònic oscillator 

( ■ I '• ') 1 l Xq 2 + q' 2 )cos^-^)-2qq' \ 

' " ' " /^7r(l-exp{-ï2(v9-^)}) eXP V 2hsin ^ ~ ^ i ' 



In the limiting case tp' — > (p one obtains the delta function S(q-q') and in case of (pf = <p+n/2 
the scalar product (q\p) if one substitutes q' = p. 

The Wigner quasiprobability contains the complete information about the state de- 
scribed by the density operator. The reconstruction of the density operator from the Wigner 
quasiprobability can be made by transition operators. As was already mentioned the Wigner 
quasiprobability has the unique property among all quasiprobabilities that the transition op- 
erators for the transition from the density operator to the quasiprobability by forming the 
trace and from the quasiprobability back to the density operator by a phase-space inte- 
gral are the same. The reconstruction formula of the density operator from the Wigner 
quasiprobability by means of these transition operators T (q,p) in real or T Q (a, a*) in com- 
plex representation has the following form 

g = 21x71 J dq A dpW(q,p)T (q,p) = tc J -da A da* W(a, a*)T (a, a*), 
W(q,p) = (gT (q,p)), W(a,a*) = (gT (a,a*)). (4.8) 

The transition operator for the Wigner quasiprobability is essentially the displaced parity 
operator [64-66,38]. There exist many representations for this transition operator. We 
use for our purpose the representation explicitly contained in Eq.(3.1). The reconstruction 
formula (4.8) can be written then in real representation in the following form 

q = 2Tití J dq AdpW(q,p)exp(^-Q-^-- P-^j 5(q)5(p) 
= 2Tin J dqAdp5{q)5{p)expÍQ^ + P-^-\w(q J p) 



^of^o m [dq k dp l 



(4.9) 

(<7=0,p=0) 



where S{. . .} means symmetrical ordering of the content in braces. This formula shows 
first of all that the Wigner quasiprobability as all other quasiprobabilities too contains a 
great redundancy and that for the reconstruction of the density operator it is already suffi- 
cient to know the Wigner quasiprobability W(q,p) in an arbitrarily small neighbourhood of 
(q = 0,p = 0). One could introduce in Eq.(4.9) the representation of the Wigner quasiprob- 
ability W(q,p) by its Radon transform W(u, v;q) explicitly given in Eq.(2.13) and has a 
representation of the density operator by the rotated quadrature components. Such a form, 
however, is not very convenient, for example, for the calculation of the matrix elements of 
the density operator in the Fock-state representation. It is better for some purposes to use 
the Fourier transform W(u, v) of the Wigner quasiprobability as an intermediate step. 
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The introduction of the Fourier transform of the Wigner quasiprobability according to 
Eq.(2.7) into Eq.(4.9) yields 

q = 2hn J dq A dp5(q)5(p) exp (Qt^ + J^y j du A dvW(u,v) exp {i(uq + vp)} 

= — íduAdvexp{i(Qu + Pv)}W(u,v). (4.10) 
2ir J 



Now, by expressing the Fourier transform of the Wigner quasiprobability by the Radon 
transform of the Wigner quasiprobability one finds ( / is the identity operator ) 

^ r + OO r 

g — — / dc / du A dvexp{-i(Ic - Qu - Pv)}W(u, v; c). (4.11) 

2,71 J-oo J 

After the substitutions 

u = r cos p v = rsm<f, c = rq, (4-12) 

one arrives at ( note the general identity J +o ° dr r J 2lT f(r, ip) = f^™ dr\r\ J W f(r, </?) for unique 
functions of r and ip over the plane ) 

fï /"7T f + OO f + OO ^ 

q — — / d(p dq dr\r\ exp{— ir(Iq — Q cos <p — P sin p)}W(cosp, sin ep; q) 

2tï JO J-oo J-oo 

1 f+°° v d ( 1 \ 

= *7o ^loo " g ^ (C0S ^ SÍn ^ g) ^r7^g^)J 

= -;Jt iq { V 7^W)) ^n^-P,^), (4.13) 

where Q(<p) is the rotated operator Q according to Eq.(3.10). Recali that V(l/x) = 71(1 /x) 
means the principal value of 1/x and its first derivative — 7Z(l/x 2 ) the canonical regulariza- 
tion of —1/x 2 as explained in section 2, Eq.(2.15). However, the operator character of the 
singularities in this formula makes it more difficult to intèrpret them. A possible way to do 
this is to make a Taylor series expansion in powers of Q(<p) in the following way 

1 /"7T r + OO l 

Q = / d p / dgftW(cosp,siny>;g)7£— — 2 

n Jo J-oo (Iq — Q((p)) 

= ~J2( k + 1 )- d<P[Q(<p)) / dg^W(cosy,siii^;g)7e-2^ 

= -f;(fc+i)i rdp(Q(p)) k 

k=0 71 Jo 

r+oo Tl ( ^ v 

x / cíg l^(cosy2, siny2; g) + (— l) A W(cosy9, siny9; — g) 

fe „z r)'W7 ~1 

-E(l + (-l) fe+ |^r(cos^sin^;0)|. (4.14) 
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In the last part of this representation we have explicitly written down the meaning of the 
canonical regularization of powers l/q 2+k ( see also Jïí| ). The integration over q goes here 
only from zero to plus infinity. We wrote this equation in detail because the content of 
formulae such as (4.13) was not represented with clarity in the literature about quantum 
tomography. However, it seems that Eqs.(4.14) does not provide a very convenient approach 
and mostly it is better to make the integration with respect to r in Eq.(4.13) not before 
calculating the matrix elements of the density operator in Fock-state representation or in 
other representation. Nevertheless, it cannot be excluded that sometimes some initial terms 
of the Taylor series expansion in (4.14) could give already a good approximation that must 
be investigated but the same approximations can be obviously obtained by Taylor series 
expansion of the exponentials in (4.13) where the integration over r is accomplished in a 
later stage of the calculations. 

A more convenient form of the reconstruction formula of the density operator from its 
Radon transform for many purposes can be obtained when the integration with respect to r 
is not accomplished before making the transition to normal ordering. Our method is similar 
to the method used in |TJ], Fw]] . In such a way by transition to normal ordering one obtains 



i /""" r+oo u d ( fi /"+ 00 T 

Q = — / dip dq W ( cos ip,smp;q) — [i— dr — exp i—irílq — QUp)) 

7T Jo J-oo ' Oq \ 2 J-oo l r l , 

dq iy (cos <p, sin <p; q) 





exp \ir\—a' exp 




(4.15) 



and by the substitution of r — > id/dq in parts of the integral over r 

j /"7T r+oo ^ 

q = — dtp dqW (cos tp, sin p; q) 

7T JO J-oo 



:v^||^)exp (-^yfl) 
\fh f +oc , r ( . h 2 



xi / dr — - exp — iqr r >. (4.16) 

2 J-oo \r\ \ 4 J J 

We obtained a representation where a normally ordered operator part is separated from an 
integral onto which it acts. Now, we make the variable substitution 

y/K' dx " /Tl dq ^ y/K' { J 



The integral over r in Eq.(4.16) can be transformed and solved in the following way 

y/fï r + oo y / fi \ f \ f x / \ 

i—^- J dr — exp í — ivTixr — — r 2 1 = 2 exp y— x 2 J J dt exp [t 2 J = 2F(x) 
= V2D (V2x) % - {D^ (iV2x) - (-iV2x) } = h (x)g (x), (4.18) 
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with the abbreviations 



ho(x) = — exp -— = —D Q (V2x) , / dx {h (x)f 

7T4 V Z J 7f 4 v ' J — oo 



# (X> = 2ttï exp (^-yj f* dt ex P ( í2 ) = v^tt^ \D_ X (iV2x) - D_i (-iv 7 ^)} , 

W(/ioOr),í/o(z)) = 2, = f(x)gW(x) - f {1 \x)g(x), (4.19) 

where W (f(x), g(x)) denotes the Wronskian of two functions f(x) and g(x) and D u (z) the 
functions of the parabòlic cylinder, in particular 

D_, (±iV2x) = y|exp (-Çj (l=F *(«)), *(*«0 = jf ds exp V) , (4.20) 

with $(z) as the error function. The error function of imaginary argument has a close 
relation to the Dawson integral denoted by F(z) in (4.18). The function ho (q/VÏÏ) /h 1 ^ 
is the normalized solution of the wave equation for the harmònic oscillator in the ground 
state that means the eigenstate of the number operator to the eigenvalue n = in "position" 
representation. It possesses even parity. The function g (q/ y/hj is under all nonnormalizable 
eigenfunctions of the number operator to the eigenvalue n = the eigenfunction with odd 
parity. 

The normally ordered operator part in Eq.(4.16) is convenient for the calculation of 
matrix elements in the Fock-state representation and in the coherent-state representation. 
We begin with the simpler last case. Due to the relation 

M exp (-^-) exp (-^-j (%))' \a) = exp ( j , (4.21) 

it acts onto the integral in (4.16) as a displacement operator of the argument and one obtains 
from (4.16) immediately the following reconstruction of the coherent-state quasiprobability 
from the Radon transform of the Wigner quasiprobability 

n , *x _ (a\g\a) 
Q{a,a ) = 



7T 



l r-TT r+OO w 2y?l O 

= — dtp dq W(cos(p,smLp; q) —F 

7T JO J-oo ix Oq 



dq v/2 J ' 

(4.22) 

Herein, F(z) denotes again the Dawson integral which definition can be taken from Eq.(4.18). 

We now consider the reconstruction of the Fock-state matrix elements. First of all, we 
find for the Fock-state matrix elements of the operator part in (4.16) 
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exp {i(m — n)ip} ^ T ^^ m\n\ f 1 d \ m+n 2j 

h j!(m- j)!(n- j)! ^"T^J 
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n' / I <9 \ / 1 d 

= «p {i(m - nMl /-j(-__) (4.23) 

where L™ _n (;a:) denote the associated Laguerre polynomials. By forming the Fock-state 
matrix elements in (4.16) and by inserting then (4.23) one obtains the following structure of 
the reconstruction formula 

1 /•■* r+°° ^ / q \ 

(m\g\n) = - dip J dq W(cosíp, sin tp; q) exp {i(m - n)íp}F m<n l~ ^1 > (4.24) 

with functions -F m , n (a;) defined in the following way 

d í 1 {m '" } m\n\ ( 1 5\ m+n - 2 ^' 

^.ní») = -K- { i== Yl 17 ^77 1 — miT h o{x)g {x) 

dx [ VmR f^ j\{m - j)\{n - j)\ \ y/2dxj 

1 m\n\ ( 1 d \ m+n ^ 

E ^ ^ „-m ~ 7^1 F o,o(x) = F„, m (*).(4.25) 



The functions exp {i(m — n)(p}F mtn (x) are called "pattern functions" for the reconstruction 
of the Fock-state matrix elements [12-14]. However, contrary to the cited papers we have 
defined them with a factor 7r larger because then the integral over the angle (p in (4.24) with 
1/tt in front can be considered as an avering over the angle. The angle-dependent part is 
explicitly splitted in (4.24) from the more complicated position-dependent part. The pattern 
functions F m>n (x) are symmetrical in the indices by definition (4.25). 

Let us consider the explicit representation of the pattern functions by functions of the 
parabòlic cylinder D v {z). Due to the relations 

(±75 ( X - S))" ± ^ = J>-n+l) J ™( ± < ^*)· (426) 

following from bàsic definitions of these functions one obtains by accomplishing the differ- 
entiations in (4.25) and taking into account (4.18) and (4.19) 



■m+n 



2 / f^o j ] -(m - j)\(n - j)\ 

{D_ (m+n+2 _ 2j) (iV2x) + (-l) m+íl LL (m+n+2 _ 2i) ( - iV2x)} . (4.27) 



This is near to the representation derived in [H] but the calculations are made there more 
generally for the reconstruction of the matrix elements in the Fock-state basis from the 
Radon transform of the s-ordered quasiprobabilities. In section 6 we will obtain an alter- 
nat i ve representation of the pattern functions in form of a series over even or odd Hermite 
polynomials. 
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We now derive a specific structure of the pattern functions as derivatives of the products 
of a normalizable with a nonnormalizable wave function of the harmònic oscillator [18-21]. 
The operator in (4.23) acts onto the product h (x)go(x) as an excitation operator. The 
creation operator can be substituted in application to wave functions (q\ip) = ip(q) by 

Therefore, the higher excitations of the eigenfunctions of the number operator can be defined 
by 

(\n),\n)') = ÍL=(|0),|0)') — (h n (x),g n (x)) = -j== [x - — J (/*,(*),«>(*)) , 

/+oo 
dx h m (x)h n (x) = 6 mtn , W(h n (x),g n (x)) = 2, (4.29) 
-oo 

where \n)' is the abstract notation for the nonnormalizable eigenstates of the number operator 
with parity (— l) rí+1 . The functions h n (x) are the orthonormalized Hermite functions or with 
the substitution x = q/VK and the absent normalization factor the wave functions 

(q\n) given in (4.5). In the same way, the annihilation operator a in "position" representation 
can be represented according to 



Q + iP 
a = — = • 

However, there is a very important difference concerning the action of the operator a to the 
vacuum state |0) and to the state |0)', in "position" representation 

-L (x + ^ h (x) = 0, ±(x + ^ 9o (x) = .^exp (Ç) = g^x), 

(4.31) 

Obviously, by repeated action of the annihilation operator in position representation onto 
the state g~i(x) one can obtain nonnormalizable eigenstates g- n {x) of the number operator 
to the eigenvalues — n whereas for the normalizable eigenstates the corresponding series is 
truncated below the vacuum state h (x). We do not investigate this here more in detail but 
show that the asymmetry in the action of the annihilation operator onto h (x) and go(x) 
affects an assymmetry in the representation of the pattern functions by the series h m (x) and 

9n(x). 

We now prové that 
ò 

F m , n (x) = — {h rn (x)g n (x)} , (m < n + 1), ( F m ^(x) = F n ^ m (x) ). (4.32) 
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The proof of this formula can be made by complete induction using the recursion relations 

S m ,n+i(u) = J—r \ ~e~ l<f -j=S mtn (u) + VnS m ^ ltn (u) \ , 
\Jn + 1 [ V2 I 



Vm + 1 1 V2 



?m+l,n(") = 7 =T -e^-TS'S'm.nÍM) + V^m.n-lH , (4.33) 



following from (4.23). The relation (4.32) is obviously true for m = and arbitrary n because 



(4.34) 



9 1 / 9 \ 2 9 



•9a; v / 2"ní \ cte / dx \Jl n n\ \ dx ) dx 

On the other side, however, one finds 

dx^¥2\\ dx) { h o(x)9o(x)} = -^{h(x)g (x)-V2x}. (4.35) 

This means that (4.32) is not generally true for all (m,n = 0) without restriction. Let us 
give the proof of (4.32) by complete induction. Suppose that it is true for a certain (m,n). 
Then by using the recursion relation (4.33) one can prové without difficulties that it is also 
true for (m, n + 1). We do not explicitly write this down because we will write down the 
induction from (m,n) to (m + l,n) which is very similar but shows in addition a specific 
difficulty By using (4.33) one finds by induction from (m,n) to (m + l,n) 

^ S m+l,n (J^j {h (x)g (x)} 

9n{x)— [x - — h m (x) - h m (x) [—= [x + — g n (x) - y/ng n ^{x) \ 



dx + 1 ( \ 9x / 

= e l(m+1 " n) ^ {h m+1 (x)g n (x) - 5 nfi h^g^x)} , (4.36) 

where the relations (4.29) and (4.30) are applied. The sum term proportional to 5 n $ in 
brackets which have its origin in the relations (4.31) prevents the full proof of (4.32) for all 
(m + 1, n). However, if one raises at once n to n + 1 that means if one makes the complete 
induction from (m,n) to (m + l,n + 1) then the proof makes no difficulties because of 
<Wi,o = f° r & U nonnegative n. This means that the first relation in Eq.(4.32) is proved 
with the restriction m < n + 1 because it is right for (m = 0, n) and (m = l,n). If we do 
not consider the restriction m < n + 1 then one obtains new pattern functions by the first 
of the relations in Eq.(4.32) which do not coincide with F m ^ n (x) defined by (4.25) or (4.27) 
for m > n + 1 but can be used by full rights as equivalent pattern functions. We investigate 
this nonuniqueness of the pattern functions for the reconstruction of the Fock-state matrix 
elements furthermore in sections 6 and 7. The pattern functions given by (4.25) or (4.27) 
possess the advantage that they vanish in infinity. 
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5. Reconstruction of the normally 
ordered moments from the Radon 
transform of the Wigner 
quasiprobability 



The normally ordered moments (a* a l g) can be directly reconstructed from the Radon trans- 
form of the Wigner quasiprobability without the intermediate calculation of the matrix el- 
ements in the Fock-state representation as found by Richter [23]. A very simple derivation 
of this reconstruction is given in |25] where is also shown that the integration over the angle 
in the Radon transform made in EM is not necessary and can be substituted by certain 
summations over discrete angles that leads to more bàsic formulas. 

The starting point for the reconstruction of the normally ordered moments is the following 
representation of the density operator by the normally ordered moments derived in [63-65] 

£ = EÍ>M(« tfc ^ a kJl = { ^ ^L·ÉÈZÉ . (5.1) 

fc=0Z=0 j=0 j\yj(k - - j)\ 

It can be easily obtained from the normally ordered expansion of the bàsic Fock-state oper- 
ators |j7i)(n| in powers of the boson operators fHf . For the Radon transform of the Wigner 
quasiprobability in the form W (cos ip, siny?; q) one finds from (5.1) 

oo oo 

W(cos(f,sm(p;q) = (q; tp\ g\q; <p) = ^^(q; p\a Kl \q; p) (a) k a l g) , (5.2) 

fc=0 1=0 

and with the explicit form of the position representation of the number states by Hermite 
functions 

(q; (p\a kt i\q; (p) = ~Í=exp 

v 7ih 
l 

exp 




irh 

Here we used a known identity for finite sums over products of Hermite polynomials ( see 
Eq.10.13. (36) in [B3] ) which can be proved by complete induction. By multiplication of 
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Eq.(5.3) by H n (q/y/nj and by using the well-known completeness relations for the Hermite 
functions one finds 

/ + OO w / fl \ _ 71 T)l 



= (jV{ (<*-** + V»)"}e) = (Jï ) (Af{(«M)"} 8 ). (5.4) 

The symbol 7V{. . .} denotes normal ordering of the content in braces. The system of equa- 
tions (5.4) for the normally ordered moments (a^ k a n ~ k g) can be solved by using the properties 
of the solutions of the circle division problem. The solutions z for the division of the unit 
circle in n + 1 equal sectors ( harmònic division ) 



= z n+1 - 1 = (z- 1) z = exp(is^-) , s = 0,1,...,ti, (5.5) 

m=o ^ n + 1 



possess the property 



n n ( 27T 

E zm =E ex P — rr ) = ( n + !) ( 5 - 6 ) 



m=0 m=0 



n+1 



By using these orthogonality relations of the solutions of the circle division problem one 
finds from Eq.(5.4) 

/ tfe l \ m í n / fmr \1 1 

{a aQ) = (FhttïíS^i ,(Jfc "° Kfh+t)}^ 

x £j rfg iy {cos + ^^ T ) , sin + ; g} . 

(5.7) 

Herein, ip denotes an arbitrary initial angle. This formula is the only place we know where 
our notations of the Radon transform seem to be a little bulky. By integration over the 
arbitrary initial angle (po one obtains [23], 25 1 



fc'i' 1 f n 1 r+oo u i q \ 

(at a e > = 7^])?- I *P ex P H(* - *M "yp^ *W (cos^, sin^; g) ^ J . 

(5. 



It is interesting to mention that from a technical point of view it is mostly more simple to 
make first the integration over the variable q in this formula for simple cases of Radon trans- 
forms such as, for example, for squeezed coherent states whereas the direct integration over 
the angle of the Radon transform multiplied by phase factors is connected with considerable 
technical difliculties. 

It arises the question whether solutions of Eq.(5.4) with unequal divisions of the unit circle 
are possible or not. Such solutions are, in principle, possible if one uses n + 1 inequivalent 
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discrete angles ( pairs of angles ip and ip + n are equivalent ). However, it seems to be difficult 
to find general solutions for arbitrary n in explicit form. For low orders of n it is possible to 
solve the system of algebraic equations obtained for discrete angles. In such a way one finds 
for the first-order moments 

1 f +co j ^ exp(i(fi)W(cos(po, sin (p ; q) exp(i(p )W(cosipi, sin</?i; q) \ 



{ag) = ^7=2 L dq 



1\pÏJ-oo \ i sin(</?i - <fo) i sin(</? - <fii) j 

xH^-^j, (<f e ) = (a e )*, H 1 (x)=2x, (5.9) 

and for the second-order moments 

2 v 1 f +oc , í exp + ip 2 )) W(cos ipo, sin ip Q ] q) 



ui: 



8 J-oo [ sin(</?i — 



(a^ag) = -- / 

O J— oc 



(fio) sin(v? - ^2) 
exp (i(y? 2 + v? )) VF(cos sin ip t ; q) 

sin(v?i - ip ) sin(</? 2 - V?i) 
exp (i(ip + (£1)) W (cos y? 2 , sin y? 2 ; g) | / g 

sin(v? 2 - sin(y? - <p 2 ) J 2 \V^ y 
1 /"+ 00 , í cos (y?2 - Vi) V^(cos(/9 , sin^ ; q) 



sin(v9i - v? ) sin(y? - y? 2 ) 
| cos (y?o - y? 2 ) iy (cos yi, sin <p^ q) 

sin(v?i - ipo) sm(ifi 2 - (fií) 
| cos (y?i - v? ) W~ (cos y? 2 , sin y? 2 ; g) \ R ( _q_\ 
sin(ifi 2 - (fií) sm(ifi - ifi 2 ) J 2 \y/%)' 

(a^g) = (a 2 ^)*, # 2 (z) = 4a: 2 - 2. (5.10) 

It is easy to specify these solutions. For example, for ifi — 0, ifií — 7r/4,</?2 = vr/2 it follows 
from (5.10) 



(a 2 Q) = l dq^(l-i)W(l,0-q)+i2W (^,y^-q^ -(l + i)W(0, l;ç)}# 2 
(aW) = Í/^rfgjvF(l,0;g) + iy(0,l;g)jiJ 2 ^^, (5.11) 
and for harmònic division </? = 0, ifií — n/3, (fi 2 = 2n/3 
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(a 2 g) 



^ dqlw(l,0;q)- 



2 




2 




■) 



I 



(or a£>) 




W(l,0;g) + W 



( 




We are not sure to what extent Eqs.(5.9) and (5.10) can be explicitly generalized. The 
harmònic division has the disadvantage that one must change the division when going from 
moments of the n-th order to moments of (n + l)-th order. If one is only interested up to 
moments of the 4-th order then one can take the harmònic division in 5 angles according to 
<fo + (mir)/5 with m — 0,1, ... ,4 and can use for the lower moments subsets of these angles 
corresponding to anharmonic divisions. In photon statistics one considers only moments 
(N l g) corresponding to linear combinations of normally ordered moments (a* k a k Q) up to the 
order /. It is quite possible that some other special solutions of the systems of equations 
(5.4) corresponding to anharmonic divisions can be found. 
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6. Alternat ive method of the 
reconstruction of the matrix elements 



The matrix elements of the density operator in Fock-state representat ion can be also re- 
constructed from the normally ordered moments of the density operator. If we use the 
reconstruction of the normally ordered moments from the Radon transform that was dis- 
cussed in the last section as an intermediate step to the reconstruction of the matrix elements 
then we arrive at a new relatively simple alternative representation of the pattern functions. 
We now consider the derivation. 

From the representation of the density operator by the normally ordered moments in 
Eq.(4.1) one obtains by forming the matrix elements with Fock states 

oo oo 

(m\g\n) = ^ ^(m\a K i\n) (a ] k a l g) 

k=0 1=0 

= v v v (-^'M* ~ M fc ~ M ( a t* i g ) 

k=0 1=0 j=0 j\y/(k - - j)\ 

1 00 ( — 



y/rrúrú j=0 j! 

By inserting the reconstruction formula for the normally ordered moments from the Radon 
transform in the form averaged over the angle (5.8) one finds after changing the order of 
summation and integration 

(m\g\n) = — d<p exp {i(m — n)<p} / dq W(cosip, sin ip; q) 

TY JO J-oo 

1 ^ (m + j)!(n + j)! ( ÍV ( q\ 

1 r-TT r+oo w / g \ 

= - J Q dtp exp {i(m - n)(f} dq W(cos(p, sin <p; q) F mn I -j= I . (6.2) 

This means that we derived a representation of the pattern functions F' mn (x) by an infinite 
series over Hermite polynomials of even or odd order in increasing steps of 2 in the indices 
with some coefficients and in dependence on the even or odd order of m + n 

* ( \ - 1 ^ (m + j)\(n + j)\ f IV 

mA ' ~ ^2^m\n\^ j\(m + n + 2j)\ \ 2) W 
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9 1 1 E i^rf·. í-^^wl· («3> 



dx [ 2V2 m+n m!n! ^ j!(m + n + 2j + 1)! V 2 

The second representat ion in (6.3) by a derivative of a function with respect to q was made 
because it corresponds to the general structure of the pattern functions as derived in [18- 
20] where the content in braces can be represented as a product of the only normalizable 
eigenfunction of the number operator with a nonnormalizable eigenfunction of this operator 
or linear combinations of such products ( superpositions of the normalizable with each non- 
normalizable eigenfunction give in every case nonnormalizable eigenfunctions, however, with 
no determined parity, in general ). We wrote i^ n (x) for the pattern functions because they 
must not necessarily be identical with the pattern functions F mjJl (x) derived in Eq.(4.27). 
We now investigate this new phenomenon of the nonuniqueness of the pattern functions. 

In the last section it was found in Eq.(5.4) that the multiplication of the Radon transform 
with a Hermite polynomial H s (q/ \^ÏÏj and its integration over q leads to a linear combination 
of (s + 1) normally ordered moments of degree s multiplied by so me binomials coeficients 
and phase factors e^ 2fc_s ^ with k — 0, 1, . . . , s. In the reconstruction formula of the Fock- 
state matrix elements from the Radon transform of the Wigner quasiprobability (4.24) the 
pattern functions F m<n (x) have to be multiplied by phase factors Q l ( m - n )^ a nd integrated over 
ip from to 7r. Because of 

1 f w 

- dif exp(i2lip) = 6i t0 , (6.4) 



7T JO 

two pattern functions F mtn (x) and F^^x) connected by 



F L,n( X ) = F m,n( X ) + C k H\ m _ n \_ 2k (x) , (6.5) 

k=l 

where [u] denotes the integer part of v and where Ck are arbitrary coeficients, lead to the same 
result for the reconstruction of the Fock-state matrix elements {m\g\n) and are equivalent. 
It is not proved that (6.5) is the most general form of the nonuniqueness of the pattern 
functions but it seems so and it is the only form of the nonuniqueness which plays a role 
in the present paper. The pattern functions obtained in (4.27) become identical with the 
pattern functions in (6.3) if we only formally extend the summation to all possible negative 
summation indices j, i.e. 

F m ,n{x) = = , oTÜ í~ïï) H m+n+2j(x) 



V2 m+n m!n! J= _[^±n] + n + 2j)\ \ 2 



F m,n( X ) 



m — n 
I 2 



(k- l + {m,n})\(\m-n\ - k)\ k 



V2\ m ~ n \m\n\ fi (k - 1)\ (\m - n\ - 2k)\ 

(6.6) 
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where {m,n} = Min(m, n). All investigated special cases are in agreement with this state- 
ment but the general proof is not made up to now. In special cases one finds 



x 



Fn,n(%) F n ^ n (x), F n+ i^ n {x) F n+1 ^ n { 

F n +2,n{x) = F' n+2tn {x) - J^^^ H (x), F n+3 , n (x) = F' n+ ^ n {x) - J ^ H x {x). 

(6.7) 

Let us discuss some properties of the pattern functions F^^x). First of all, they are 
symmetric in the indices and possess the parity (— l) m+n as F mn (x) too 

KJ?) = F^ m {x), F' m>n {-x) = (-l) m+n F^ n (x). (6.8) 

The summations at the point x = can be exactly accomplished for not very large even 
differences \m — n\ whereas for all odd differences the corresponding Hermite polynomials in 
(6.3) vanish for x = 0. One obtains [ 7^ ] 

In 4- 3 

<,n(0) = (-l) n 2, F' +2 , n (0) = (~1)" +1 = , F' n+2l+l ,M = 0. (6.9) 

y/(n + 2)(n + l) 

The first màxima and minima of F' nn {x) = F n ^ n (x) for ï / come for increasing n very 
near to the vàlues ±2 but are not exactly equal to these vàlues |72[. Some similarities in the 
form of the pattern functions to the corresponding products of wave functions (q\m)(n\q) 
can be explained by the following expansion of products of Hermite polynomials which can 
be proved by complete induction ( see Eq.10.13 (37) ) in |63[ ) 




(q\m)(n\q) = . — exp — — =H„ 
^fï V n V2 m + n m\n\ 




{m,n} , , 

H m {x)H n {x) = -Tt ^77 ^H m+n _ 2j (x). (6.10) 

p ]\(m - ])\{n - ])\ 

The first polynomial terms in the expansions in (6.3) and (6.10) are both H m+n (x). Contrary 
to (6.3), the next terms in the expansion in (6.10) are proportional to Hermite polynomials 
with indices decreasing in steps of 2 but if the normalized Gaussian function exp(— x 2 ) / y/ïr is 
included into the expansion in terms of Hermite polynomials, then we get a more complicated 
formula which we do not derive here. Thus one has similarities in the representation of the 
products of wave functions and corresponding pattern functions by expansions in Hermite 
polynomials. 

The representation (6.3) of the pattern functions is appropriate for their calculation by 
a computer with good accuracy for not too large arguments ( say \x\ < 5 with ~ 50 initial 
sum terms and not too high (m,n), e.g., for m = n = one obtains by 15 initial terms an 
approximation where the deviations begin for \x\ ~ 5 ). For large vàlues of the argument 
one can use the highest power in the Hermite polynomials H n (x) that means (2x) n as an 
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approximation. It is interesting that the Hermite polynomials H n (x) can be obtained from 
these asymptotic functions (2x) n by the following convolution [[F3|, [74j] 



H n {x) = exp (2x) n = ^exp (x 2 ) * (2x) n . (6.11) 

One can prové by Taylor series expansion of the exponential function that one obtains by 
this convolution the explicit representation of the Hermite polynomials. By applying this to 
Eq.(6.3) one finds 

' 1 d 2 \ (V2x) m+n ~ (m + j)\(n + j)\ / n 2 y 



FUx)=ew {--^ (-2X»)-. (6,2) 

If we omit the convolution operator exp (— (1/A)d 2 /dx 2 ) in this relation we get an asymptotic 
representation of the pattern functions for large vàlues \x\. 

If we use the explicit representation of the Hermite polynomials in the pattern function 
Fqo( x ) = ^o,o ( x ) we arrive, after reordering of the sum terms and accomplishing one of the 
sums, at the following Taylor series representation 

d i ^ j\ ( iy 

- 2g W (2I)2 '' (613) 

This can be also represented in the forms 

d 



Fq,q(x) = — í^exp^x 2 ) £ dt exp(í 2 ) 



(V f V2 exp (-Ç ) l -{D-t (iV2x) -iVfcir)} 



&r V 2/2 



= exp (-^jj {D- 2 (iV2x) + D_ 2 ( - iV2x)} . (6.14) 

where D v (z) denotes the functions of the parabòlic cylinder. Let us show that the content in 
braces in Eq.(6.3) is proportional to the product of the normalizable eigenfunction of order 
m and the nonnormalizable eigenfunction of order n. For this purpose we make the following 
transformation 



( d\ m ( x 2 \\ { ( d \ n _ ( x 2 " 
x — — 



exp (" Y j } { [ x ~ Yx ) 2 exp (" Y ) I* du exp M 

2H m (x) (—7^ J { exp ( _a;2 )/ d«exp(w 2 )j 
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k=0 



1 « (m+j)!(n + j)! ( IV 

Jt7~~ , ^ , o.- . im ~o H m+n+2j+1 (x). (6.15) 



2 3 . p„ +1] j!(m + n + 2j + l)! V 2 

In the first step we brought the operator {x — d/dx) n from the right to the left of the function 
exp(— x 2 /2). Then we applied (6.13) and substituted the product of Hermite polynomials 
by sums over Hermite polynomials according to a known formula ( p3| , 10.13 (36) ) which 



can be proved by complete induction. After rearrangement of the double sum one sum is 
calculated by using 

™ (-iy m \( n + j + i)\ ( m + J )!( n + J )i 

f^l\(m-l)l(n + 2j + l+t)\ j\(m + n + 2j + l)V [ ' ' 

which can be proved again by complete induction ( m — > m + 1 because it is true for m = 
with arbitrary n and j ). Thus there is directly proved the connection between two essen- 
tially different representat ions of the pattern functions. We think that the representation 
of the pattern functions by infinite sums over Hermite polynomials possesses some advan- 
tages in comparison to the representation by sums over functions of the parabòlic cylinder 
of imaginary argument. 



34 



7. Differential equations for products 
of Hermite functions and 
orthogonality relations 



Initiated by Richter |I8J to explain the structure of the pattern functions as first derivatives 
of the product of the normalizable and one nonnormalizable eigenfunction of the number op- 
erator in "position" representat ion it was found [[[1], that the most appropriate approach 
to solve this problem is to consider the differential equation for products of Hermite functions 
and to establish the orthogonality relations with the help of the adjoint differential equation. 
This will be presented here. The problem solved in |Tj| was from one point of view more 



general because there were considered products of eigenfunctions of more general Hamilton 
operators for one degree of freedom consisting of the usual part from kinetic energy and of 
a part from arbitrary potential energy but from the other point of view it was more special 
because there were considered only products of eigenfunctions to the same eigenvalue. This 
last restriction is omitted in ||21|| . The most interesting case for the nondiagonal elements 
corresponding to products of eigenfunctions to different eigenvalues is, however, the case of 
the harmònic oscillator because for the most other interesting one-dimensional potentials 
defined over the whole coordinate axis R one has no degeneracy of energy differences in the 
eigenvalue spectrum. The multiplication of the Radon transform with phase factors and 
integration over the angle in case of the harmònic oscillator sorts out the Fock-state matrix 
elements corresponding to equal energy differences of the energy eigenvalue spectrum. 

The differential equations for both series of normalizable and nonnormalizable excitation 
states h n (x) and g n {x) of the harmònic oscillator in position representation are ( as before 
we use the abbreviation x = q/ \/h ) 

í d 2 2 „ A „ , , „ (x 2 l d 2 



dx 2 



x 



2n + lL·(*)=0, — í^-i_j /„(*)= (n + i) /„(*), (7.1) 



where f n {x) is written for an arbitrary of the functions h n (x) or g n (x) or their linear com- 
binations. The normalizable functions h n (x) are the Hermite functions. They possess the 
parity (— l) n and are given up to complex normalization factor in a unique way and are 
connected with the Hermite polynomials H n (z) and the functions of the parabòlic cylinder 
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D n (z) as follows 



1 í x 2 \ H n (x) 
D n (y/2x) 



7T4 vn! 



\pïn\ i 



- {D^ n {iV2x) + (-l) n LU-n( - ïv^x)} . (7.2) 



7T4 2 

They satisfy the following orthonormality and completeness relations 

/+oo °° 
dxh m (x)h n (x) = 6 mtn , K(x)h n (y) = S(x - y). (7.3) 

"°° n=0 

All other linearly independent solutions of equation (7.1) are nonnormalizable but among 
them one can choose g n (x) in a unique way as the solutions with parity (— 1)™ +1 defined 
in (4.29). The explicit form of these solutions expressed by the functions of the parabòlic 
cylinder is given by 

■i n+1 



where F(x) denotes the Dawson integral defined in (4.18) ( see also (4.20) for the connection 
to the error function of imaginary argument ). From (7.2) and (7.4) one finds that functions 
proportional to (iy/2x^j and D_i_ n (— i\f2x^j could be also chosen as two linearly 

independent solutions of the equations (7.1) but they are both nonnormalizable and do not 
possess definite parity. 

Let us now consider the differential equation for the product f m (x)g n (x) of two functions 
f m {x) and g n (x) both satisfying differential equations of the form (7.1) ( g n (x) is here not 
necessarily the special eigenfunction with parity (— l) n+1 ). Then one finds up to second- 
order derivatives 

d 

0^{fm{x)g n {x)} = f£\x)g n {x) + f m (x)g n 1) (x), 

{f m (x)g n (x)} = fiï\x)g n (x) + 2f£\x)gW(x)+f n (x)g n 2 \x) 

= 2f^{x)g n 1 \x) + 2(x 2 -{m + n + l))f m {x)g n {x). (7.5) 
In the next step by further differentiation and by applying (7.1) one finds 
í d 3 



2 (^(x 2 - (m + n + 1)) ^- + ^(x 2 - (m + n + 1))^ | f m (x)g n 



(x) 



{ dx 3 

-2(m - n) {f m {x)g£\x) - f£\x)g n (x)} , (7.6) 
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written in a symmetrical form with respect to the derivatives. It is a pure differential equation 
of third order for f m (x)g n (x) only in case of m — n. This is understandable because in this 
case we can form 3 linearly independent combinations of the bàsic solutions, for example, 
(h n (x)) 2 , h n (x)g n (x), (g n (x)) 2 if h n (x) and g n (x) are two linearly independent solutions of 
(7.1). For m ^ n one needs a further differentiation of this equation to obtain the following 
fourth-order differential equation 

{^- 2 è(( a;2 - (m + n + 1) )|: + |:( a;2 - (m + n + 1) )) +4(m - n)2 } 

*f m (x)g n (x) = 0. (7.7) 

We have here 4 linearly independent solutions of this differential equation as which can 
be chosen, for example, the products h m (x)h n (x),h m (x)g n (x), g m (x)h n (x), g m (x)g n (x). The 
operator in the differential equation (7.7) is neither selfadjoint nor anti-selfadjoint. The 
adjoint differential equation to (7.7) for functions X mjTl (x) is 

{^ _2 (^ 2 ~ (m + n + 1) ^ + ^^ 2_(m + n + 1) ^ +4(m " n)2 } 
xX m , n (x) = 0. (7.8) 

If we differentiate this equation once more then we get immediately 

{ J? - 2 í (V - (m + - U x2 - (m + - + •>)) + 4(m - n)2 } 

x^i„W=0. (7.9) 

This means that {d/dx)X m ^ n (x) satisfies the differential equation (7.7) for products of func- 
tions f m (x)g n (x) and their linear combinations. 

We now make the proper derivation of the orthogonality relations. For this purpose we 
first write down the following two differential equations according to (7.7) and (7.8) 

{ - 2 é (V - (m + " + + U x * - (m + - + '>)) + 4(m - n)2 } 

x Y m ^ n (x) = 0, 
\dx 4 

xX kjl (x) = 0, (7.10) 



2^-(k + l + l))^ + ^( X 2 -(k + l + l))^+A(k-l) 



2 



where Y miTl (x) is written for arbitrary products f m (x)g n (x) of solutions of (7.1) or their 
linear combinations. If we multiply the first equation with Xk,i{x) and the second equa- 
tion with Y mjTl (x) and subtract the obtained equations then we find the following possible 
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representation of the resulting equation 

= £.{x kJl (x)Yj®(x) - X<8(x)Yff n (x) + X$(x)YÍl(x) - X${x)Y m , n {x) 

-2{2x 2 -( m + n + k + l + 2)){x Kl {x)Y^{x) - X ( Q{x)Y m ^x)} - AxX Kl {x)Y m , n {x)^ 
+2{m + n -k-l) {x hJl (x)Y<a(x) + X§{x)Y m , n {x)} 

+4 ((m - nf — (k - if) X Kl {x)Y rn , n {x). (7.11) 

By integration of this equation over the whole coordinate axis R under the assumption that 
the products X^}{x)Y^ n {x) vanish for x — > ±oo one finds 



( m + n - k - l) J 2 dx {x kíl (x)YÍl(x) + XÍ 2 ](x)Y m , n (x)} 
+2 {(m - nf — (k — lf} |_ + °° dx X kil (x)Y m>n (x) = 0. (7.12) 
The first of the integrals in (7.12) can be transformed by partial integration as follows 

F°dx {x kil (x)YÍl(x) + X^(x)Y m , n (x)} 

= 2 J^dx (J-X^ix^j Y m , n (x). (7.13) 

From (7.13) and (7.12) one can derive different orthogonality relations. Consider the case 

m — n = k — l = —j, — > l = k + j, n = m + j, m + n — k — l = 2{m — k). (7.14) 
Since the integral becomes nonvanishing only for m — k = one obtains 

/oo ^ (Í X ^ k+j ( X ^) Ym > m +i( X "> = ^j-oo ^ (y^ X S m +^ X ^j Y rn,rn +j (x). (7.15) 

In particular, by choosing 

x k]( x ) = h k (x)gi(x), Y mjn (x) = h m (x)h n (x), (7.16) 

where gi(x) denotes the nonnormalizable solution (7.4) of Eq.(7.1) with parity (— 1)™ +1 one 
finds from (7.15) with the help of (7.13) 

r+oo ( Q \ 

J dx \—h k (x)g k+j (x)\h m (x)h m+j (x) 
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5 i / 



+00 r / ^ \ 

^ I 1 Q^ h ^( x )9m+j(x) I h m (x)h m+j (x) 



3G 



-/i m (x)5i m+j (x) ^^/i m (x)/i m+j (x)j I 

X r+00 
= 5jfe, m -WX^m+j(z),0 m +j(2;)) / dxh m (x)h m (x) 

where the special value of the Wronskian given in (4.29) was used. The special case j — of 
the orthogonality relations (7.17) can be obtained from the third-order differential equation 
(7.6) for m = n in an easier way but for shortness we did not separately give its derivation. 

Using the representation of h m (x) by Hermite functions given in (7.2) and the representa- 
tion of g n {x) by functions of the parabòlic cylinder given in (7.4) one obtains by differentiation 
of their product the following representation of the pattern functions (4.25) 

d 



= exp f^-Çj ^J^r{(n + l)H m (x) (D. 2 ^(iy/2x) + (-l) n D_ 2 _ n (-iV2x)) 

+iV2mH m _ 1 {x) (£)_!_„ (iV2x) - {-l) n D^[ - iV2x)) j. (7.18) 

This form of the pattern functions is not symmetric with respect to permutation of the indices 
(m, n). It is identical with the pattern functions F m>n (x) in Eq.(4.25) only for m < n + 1 as 
discussed in section 4 but can be taken as an equivalent pattern function for arbitrary (m, n) 
due to their nonuniqueness as discussed in section 6. Furthermore, it is clear that, instead 
of (7.16), one can take for X^j(x) the combination with exchange of the normalizable and 
nonnormalizable solutions of the wave equation, i.e. 

x k]( x ) = 9k(x)hi(x), Y mjn (x) = h m (x)h n (x), (7.19) 

or linear combinations of them. In case of (7.19) we find as the equivalent pattern functions 

= ^{9m(x)h n (x)} = FlJx). (7.20) 

These pattern functions are identical with the pattern functions F m ^ n (x) considered in 
Eqs.(4.25) and (4.27) in the cases 

F^Jx) = F m>n {x), m<n+l, F^ n (x) = F m>n {x), n < m + 1, (7.21) 

and are connected with the pattern functions (x) considered in the last section by 

K,n( X ) = \ (K,n( X ) + Kmi*)) ■ (7.22) 
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Nevertheless, in all other cases of nonidentity one has completely equivalent pattern functions 
which can be used by same right as the pattern functions F m ^ n {x). As already mentioned, 
the pattern functions F mjn (x) possess the advantage that they vanish at infinity. In the 
special cases m = n and m = n±lwe obtained the same pattern functions by all considered 
methods. 

The nonuniqueness of the pattern functions may have as the cause a certain redundancy 
in the information contained in the Wigner quasiprobability and in its Radon transform 
but it seems to be dimcult to eliminate this redundancy. The pattern functions themselves 
play an auxiliary role and only the integrals (4.24) over the pattern functions leading to the 
Fock-state matrix elements possess an invariant meaning. The representations of the pattern 
functions in the form (7.18) or (7.20) is useful for computer calculations. In comparison to 
(4.27) they possess the advantage that one has to calculate here only a sum over the real or 
imaginary parts of two functions of the parabòlic cylinder of imaginary argument instead of 
maximally n + 1 such functions there. 
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8. Conclusion 



We introduced in the present paper the Radon transform of the Wigner quasiprobability 
in its more general canonical representation in comparison to the usual and considered the 
reconstruction of the density operator via the Fock-state matrix elements and via the nor- 
mally ordered moments. The inverse two-dimensional Radon transformation was considered 
in detail because some moments seem to have been clarified in the literature about quantum 
tomography. The transformation properties of the Wigner quasiprobability and its Radon 
and Fourier transforms with regard to displacement and squeezing of an initial state are 
considered with many technical details. This has a more practical aspect, for example, for 
the calculation of the Radon transform of squeezed coherent states and squeezed-state exci- 
tations which becomes much easier by application of the derived formulae. The calculation 
of the pattern functions for the reconstruction of the Fock-state matrix elements of the 
density operator via the normally ordered moments leads to a new representation of these 
functions by convergent series over Hermite polynomials of even or odd order which provides 
an alternative for calculation and plotting of these functions. 

We did not consider in the present paper, for example, the influence of imperfect measure- 
ments which lead to some smoothing of the Radon transform with the problem to reconstruct 
the density operator from these smoothed Radon transforms, a problem, now intensively dis- 
cussed in the literature. Furthermore, we did not extend the considerations to multi-mode 
cases but many features can be probably translated to these more general cases in a sim- 
ple way We hope that our considerations are useful for the clarification of some principal 
problems connected with two-dimensional Radon transforms and its application in quantum 
òptics. 
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